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ABSTRACT 

The ’’reverse-dynamo” mechanism — the amplification/generation of fast 
plasma flows by micro scale (tnrbnlent) magnetic fields via magneto-flnid con- 
pling is recognized and explored. It is shown that macroscopic magnetic fields and 
flows are generated simnltaneously and proportionately from microscopic fields 
and flows. The stronger the micro-scale driver, the stronger are the macro-scale 
prodncts. Stellar and astrophysical applications are snggested. 
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The generation of macroscopic magnetic fields (primarily from microscopic velocity 
fields) defines the standard ’’dynamo” mechanism. The dynamo action seems to be a very 
pervasive phenomenon; in fusion devices as well as in astrophysics (stellar atmosphere, MHD 
jets) one sees the emergence of macro-scale magnetic helds from an initially turbulent sys¬ 
tem. The relaxation observed in the Reverse Field pinches is a vivid illustration of the 
dynamo in action. Search for interactions that may result in efficient dynamo action is one 
of the most flourishing helds in plasma astrophysics. The myriad phenomena taking place 
in the stellar atmospheres (heating of the corona, the stellar wind etc.) could hardly be 
understood without knowing the origin and nature of the magnetic held structures weaving 
the corona. 

The conventional dynamo theories concentrate on the generation of macroscopic mag¬ 
netic helds in charged huids. With time the dynamo theories have invoked more and more 
sophisticated physics models — from the kinematic to the magneto hydrodynamic (MHD) 
to, more recently, the Hall MHD (HMHD) dynamo. In the latter theories the velocity held 
is not specihed externally (as it is in the kinematic case) but evolves in interaction with 
the magnetic held. Naturally both MHD and HMHD ’’dynamo” theories encompass, in 
reality, the simultaneous evolution of the magnetic and the velocity helds. If the short- 
scale turbulence can generate long-scale magnetic helds, then under appropriate conditions 
the turbulence could also generate macroscopic plasma hows. In this context, a quota¬ 
tion from a recent study is rather pertinent: the structures/magnetic elements produced 
by the turbulent amplihcation are destroyed/dissipated even before they are formed com¬ 
pletely (Bellot Rubio et al. 2001; Socas-Navarro & Manso Sainz 2005; Blackman 2005) cre¬ 
ating signihcant hows or leading to the heating. 

If the process of conversion of micro-scale kinetic energy to macro-scale magnetic energy 
is termed ’’dynamo” (D) then the mirror image process of the conversion of micro-scale 
magnetic energy to macro-scale kinetic energy could be called ’’reverse dynamo” (RD). It is 
convenient to somewhat extend the dehnitions - the D (RD) process connotes the generation 
of the macroscopic magnetic held (how) independent of the mix of the microscopic energy 
(magnetic and kinetic). 

Within the framework of a simple HMHD system, we demonstrate in this paper that the 
Dynamo and the Reverse Dynamo processes operate simultaneously — whenever a macro¬ 
scopic magnetic held is generated there is a concomitant generation of a macroscopic plasma 
how. Whether the macroscopic how is weak (sub-Alfvenic) or strong (super-Alfvenic) with 
respect to the macroscopic held will depend on the composition of the turbulent energy. 
We shall derive the relationships between the generated helds and the hows and discuss the 
conditions under which one or the other process is dominant. In Sec.l we display an ana- 
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lytical calculation based on the conversion of micro scale magnetic and kinetic energy into 
macroscopic fields and flows. In particular, we dwell on the reverse dynamo mechanism: the 
permanent dynamical feeding of the flow kinetic energy through an interaction of the micro¬ 
scopic magnetic held structures with weak hows (seed kinetic energy). In Sec.2 we illustrate 
that the theoretically derived processes do indeed take place by presenting simulation results 
from a general two huid code that includes dissipation. 


1. Theoretical Model Analysis 


The physical model exploited for how generation/acceleration is simplihed HMHD - a 
minimal model that entertains two interacting scales that can be quite disparate; the macro¬ 
scopic scale of the system is generally much larger than the ion skin depth, the intrinsic micro 
scale of HMHD at which ion kinetic inertia ehects become important (Mahajan & Yoshida 1998; 
Mahajan et al. 2001; Ohsaki et al. 2002; Yoshida et al. 2004). In HMHD the ion (n) and 
electron = (v — j/en)) how velocities are diherent even in the limit of zero electron 
inertia. In its dimensionless form, HMHD comprises of 


db 

'm 

dv 

'm 


V X 


V — ao'V X b] X b 


n X (V X n) -I- (V X 6) X 6 — V 



( 1 ) 

( 2 ) 


with the standard normalizations: the density n to no , the magnetic held to the some 
measure of the ambient held Bq and velocities to the Alfven velocity Vao- We assign equal 
temperatures to the electron and the protons so that the kinetic pressure p is given by: 
p = Pi + p^ cx 2nT, T = Ti cx T^. We note that the Hall current contributions become 
signihcant when the dimensionless Hall coefficient ao = Xio/Ro {Rq - the characteristic scale 
length of a system and Ajo = c/uio is the collisionless skin depth) satishes the condition: 
ao > rj, where rj is the inverse Lundquist number for the plasma. For a typical solar 
plasma, in the corona, the chromosphere and the transition region (TR), this condition 
is easily satished (ao is in the range 10“^° — 10“^ for densities within (10^^ — 10®) 
and p = /{AtiVaqRqO^) ~ 10“^"^, where Rq is solar radius, a is the plasma conductivity). 

In such circumstances, the Hall currents modifying the dynamics of the microscopic flows 
and helds could have a profound impact on the generation of macroscopic magnetic helds 
(Mininni et al. 2003) and fast flows (Mahajan et al. 2002; Mahajan et al. 2005). 

In the following analysis ao will be absorbed by choosing the normalizing length scale 
to be Ajo- Let us now assume that our total helds are composed of some ambient seed helds 
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and fluctuations about them, 


b = H + bo + b, V = U + Vo + V 


( 3 ) 


where bo, Vq are the equilibrium flelds and H, U and b, v are, respectively, the macroscopic 
and microscopic fluctuations. 

Notice that our ambient flelds are allowed to have a component at a microscopic scale. 
For analytical work, we choose for the ambient flelds a special class of equilibrium so¬ 
lutions to Eqs. (1-2). These solutions, also known as the Double Beltrami (DB) pair 
(Mahajan & Yoshida 1998), come into existence because of the interaction of flows and flelds; 
the Hall term is essential for their formation. The DB configurations are known to be robust 
and accessible, through a variational principle, for a variety of conditions including inho¬ 
mogeneous densities. Non constant density cases do display many interesting phenomena 
(Mahajan et al. 2002; Mahajan et al. 2005), but the dynamo and reverse dynamo actions 
can be very adequately described by the analytically tractable constant density system. We 
shall, therefore, choose the following DB pair (obeying the concomitant Bernoulli condition 
V(po + 'i^o^/2) = const (Ohsaki et al. 2001; Ohsaki et al. 2002)) 


bo 

-h V X ho = '^0, bo + V X Vo = dvo, 

a 


( 4 ) 


as a representative ambient state. The general solution is expressible in terms of the single 
Beltrami flelds G± that satisfies V x G{X) = XG{X): 


bo = C'+G+(A+) + C'_G_(A_), (5) 

Vo = (a-i + A+)C'+G+(A+) + (a-i + A_)C'_G_(A_). (6) 


Here C± are the arbitrary constants and the parameters a and d are set by the invariants 
of the equilibrium system; the magnetic helicity hio = f (Ag-ho) d^x and the generalized he- 
licity ^20 = /(Ag-l-t^o)■ VX (Ao-|-'yo)<^^ 2 ; (Mahajan & Yoshida 1998; Mahajan et a/. 2001); 
here Aq is the vector potential of the ambient held. The inverse scale lengths A+ and A_ 
are fully determined in terms of a and d: X± = \[{d — a“^) ± \/{d + a~^Y — 4 ]. As the DB 
parameters a and d vary, X± can range from real to complex values of arbitrary magnitude^. 

Our primary interest is the creation of macro flelds from the ambient micro flelds. Some 
what later we will assume, for simplicity, that our zeroth order flelds are wholly at the 
microscopic scale. This allows us to create a hierarchy in the micro flelds, the ambient flelds 


^ In the analysis below we will use A for the micro-scale and fj, for the macro-scale structures. 
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are much greater than the fluctuations at the same scale (|h| |6o|, |u| l-UoD- Following 

(Mininni et al. 2003), we may derive the following evolution equations: 

dtU = Ux{VxU) + VxHxH 

+ (^Vq X (V X v) + V X (V X Vq) + (V X bo) X b + (V X b) X bo 

( 

- (V(uo ■ t^)) - V p + — 


dv 

'm 

db 

m 

dH 


= -{U-V)vo + {H-V)bo, 

= (if-V)Ueo-(C/-V)ho, 

= V X /[Fe X ho] + "yeo X + V X [{U — V X H) X H], 


(7) 

( 8 ) 
(9) 

( 10 ) 


where the brackets < .. > denote the spatial averages and Ueo = Vq — 'V x bo. This set of 
equations can be regarded as a closure model of the Hall-MHD equations, which are now 
general in two respects: 1) it is a closure of the full set of equations, since the feedback of 
the micro-scale is consistently included in the evolution of both if and U ; 2) the role of the 
Hall current (especially in the dynamics of the micro-scale) is also properly accounted for 
(see (Mininni et al. 2003; Mininni et al. 2005) for details). 

We now choose the constants a and d so that the two Beltrami scales become vastly 
separated (since these constants reflect the values of the invariant helicities, it is through a 
and d that the helicities control the final results). In the astrophysically relevant regime of 
disparate scales (the size of the structure is much greater than the ion skin depth), we shall 
deal with two extreme cases : (i) a ~ d 3> 1 , (a — d)/ad 1 (A ~ d, p ~ (a — d)/ad 
), and (ii) a ~ d 1 , (a — d)/ad S> 1 (A ~ (a — a“^), fi {d — a) ). At this time, 
we would like to draw the reader’s attention to the origin of scale separation in the original 
equilibrium system - it is the Hall term that imposes the micro scale (ion skin-depth) on 
the macroscopic MHD equilibrium. 

Consistent with the main objectives of this paper, we will now assume that the origi¬ 
nal equilibrium is predominantly micro-scale (condition applicable for many astrophysical 
systems), i.e, the basic reservoir from which we will generate macro scale fields is, indeed, 
at a totally different scale. Neglecting the macro scale component altogether, the assumed 
equilibria becomes simpler with the velocity and magnetic fields linearly related as 


-^0 


— ho (a + n 


( 11 ) 


leading to 
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VeO = 

Vq — 'V X bo = bo a ^ 

(12) 

b = 

{a-^H - U) ■ Vho 

(13) 

V = 

[H -{X + a-^) U) ■ Vho- 

(14) 


Notice the preponderance of nonlinear terms in the evolution equations for U and H . One 
would expect that these terms will certainly play a very important part in the eventual 
saturation of the macroscopic helds, but in the early acceleration stage when the ambient 
short scale energy is much greater than the newly created macroscopic energy, these terms 
will not be signihcant. Deferring the fully nonlinear to a later stage, we shall limit ourselves 
to a ’’linear” treatment here. Neglecting the nonlinear terms and manipulating the system of 
equations, we readily derive (after ’’solving” for and eliminating the short scale fluctuating 
helds) 


H ~ (i-^-l^{Vx{H-V)boxbo), (15) 

\ a J 

U ~ ((A + a-i) (Ah - V X h) - (A + a-i) V(ho-h) X 6o> 

- {|(A6-V X h) X ho) • (16) 

where the spatial averages are yet to be performed. We use the standard isotropic ABC 
solution of the single Beltrami system. 


h ^0 


U ^ 

boz — 


Vs 


[sin Xy + cos A^:] , 
[sin A^ + cos Ax], 
[sin Ax + cos Xy]. 


( 17 ) 


to compute the spatial averages. After some tedious but straightforward algebra, we arrive 
at the hnal acceleration equations 


C/ 


V X 

2 3 


(AD)Al)c/-A/f 


(18) 


H 




V X if. 


(19) 


where Bq measures the ambient micro scale magnetic energy (also the kinetic energy because 
of (11)). The coefficients in these equations are determined by a and d (A = A(a,d)). 



7 


We see that, to leading order, H evolves independently of U but the reverse is not true: 
the evolution of U does require knowledge of H. 

In the dynamo context, the Hall-currents in the micro-scale are known to modify the a 
coefficient so that it survives the standard cancellation of the kinetic and magnetic contribu¬ 
tions for Alfvenic perturbations (Mininni et al. 2002). It is also known that, depending on 
the state of the system, the Hall effect (by replacing the bulk kinetic helicity by the electron 
flow helicity) can cause large enhancement or suppression of the dynamo action as compared 
to the standard MHD (Mininni et al. 2005). 

Writing (18) and (19) as 

H = -r{VxH), tj = V x[sU -qH], (20) 

where 

r = - Xa~^ - a~‘^) , s = A ^[(A-f a“^)^ - 1] , g = A^ ^ , (21) 

0 0 0 

and fourier analyzing, one obtains 

=-ir{kx H) , = ikx{sU-qH). (22) 

yielding the growth rate, 

= —\r\ (k) , (23) 

at which H and U increase. The growing macro helds are related to one another by 

U = H. (24) 

s + r 

We shall now show how a choice of a and d fixes the relative amounts of microscopic 
energy in the ambient helds and consequently in the nascent macroscopic helds U or H. We 
persist with our two extreme cases: 

(i) For a ~ d 3> 1 , the inverse micro scale A ~ a ^ 1 implying ~ a bo ^ he, the 
ambient micro-scales helds are primarily kinetic. These type of conditions may be met in stel¬ 
lar photospheres, where the turbulent velocity held at some stage can be dominant although 
some bo is present as well. For these parameters, it can be easily seen that the generated 
macro-helds have precisely the opposite ordering, U ~ a~^H H. This is an example of 
the straight dynamo mechanism. Micro scale helds with kinetic dominance create, preferen¬ 
tially, macro scale helds that are magnetically dominant — super-Alfvenic ’’turbulent hows” 
lead to steady hows that are equally sub-Alfvenic (remember we are using Alfvenic units). 



It is extremely important, however, to emphasize that the dynamo effect (dominant in this 
regime) must always be accompanied by the generation of macro-scale plasma flows. This 
realization can have serious consequences for dehning the initial setup for the later dynamics 
in the stellar atmosphere. The presence of an initial macro-scale velocity held during the hux 
emergence processes is, for instance, always guaranteed by the mechanism exposed above. 
The implication is that all models of chromosphere heating / particle acceleration should 
take into account the existence of macro-scale primary plasma hows (even weak) and their 
self-consistent coupling (see (Mahajan et al. 2001; Ohsaki et al. 2002; Mahajan et al. 2005) 
and references therein). 

(ii) For a ~ d 1 the inverse micro scale A ~ a — a~^ 3> 1. Consequently -no ~ ^ 

6o, and the ambient energy is mostly magnetic. These conditions might pertain in certain 
domains in the photospheres or chromospheres, where the turbulent velocity held may exist, 
but the turbulent magnetic held is the dominant component. This micro-scale magnetically 
dominant initial system creates macro-scale helds U ~ a~^H 3> H that are kinetically 
abundant. The situation has fully reversed from the one discussed in the previous example — 
starting from a strongly sub-Alfvenic turbulent how, the system generates a strongly super- 
Alfvenic macro-scale how; this mode of conversion could be called the ’’reverse dynamo” 
mechanism. In the region of a given astrophysical system where the huctuating/turbulent 
magnetic held is initially dominant, the magneto-huid coupling induces efficient/signihcant 
acceleration and part of the magnetic energy will be transferred to steady plasma hows. The 
eventual product of the ’’reverse dynamo” mechanism is a steady super-Alfvenic how — a 
macro how accompanied by a weak magnetic held (compare with (Blackman & Field 2004) 
for a magnetically driven dynamo. In this study magnetic held growth on much larger 
scales, and signihcant velocity huctuations with hnite volume averaged kinetic helicity are 
found). It is tempting to stipulate that ’’reverse dynamo” may be the explanation for the 
observations that fast hows are generally found in weak held regions of the solar atmosphere 
(Woo & Habbal & Feldman 2004). 

This simple analysis has led to, what we believe, are several far-reaching results: (1) the 
dynamo and ’’reverse dynamo” mechanisms have the same origin - they are manifestation of 
the magneto-huid coupling; (2) The proportionality of U and H implies that they must be 
present simultaneously, and the greater the macro-scale magnetic held (generated locally), 
the greater the macro-scale velocity held (generated locally); (3) the growth rate of the 
macro-scale helds is dehned by DB parameters (hence, by the ambient magnetic and gener¬ 
alized helicities) and scales directly with the ambient turbulent energy 6o (t^o)- Thus, the 
larger the initial turbulent (microscopic) magnetic energy, the stronger the acceleration of 
the how. We believe that these novel results will surely help in advancing our understanding 
of the evolution of large-scale magnetic helds and their opening up with respect to the fast 
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particle escape from the stellar coronae. This effect may also have important impact on the 
dynamical and continuous kinetic energy supply of plasma flows observed in various astro- 
physical systems. We would add here that in this study both the initial and final states have 
finite heliciies (magnetic and kinetic). The helicity densities are dynamical parameters that 
evolve self-consistently during the process of flow generation. It is also important to notice 
that the end product of the reverse dynamo action is a macroscopic flow (produced from a 
microscopic helical magnetic field) while for ’’inverse dynamo” (Blackman & Field 2004) it 
is still the macroscopic magnetic field but produced from a velocity field with helicity. 

We end the analytical section by a remark on the nonlinear terms in Eqs. (7,10) that 
do not appear later. It is amazing that the linear solution given in Eqs. (22-24) makes the 
nonlinear terms strictly zero. Thus the solution discussed in the last section is an exact (a 
special class) solution of the nonlinear system and thus remains valid even as U and H grow 
to larger amplitudes. This interesting but peculiar property that a basically linear solution 
solves the nonlinear problem pertains to both MHD and HMHD. In MHD, for example, 
it manifests itself as Walen’s nonlinear Alfven wave (Walen 1944a; Walen 1944b) while in 
HMHD it is revealed through the recently discovered solution of (Mahajan & Krishan 2005). 


2. A Simulation Example 

In order to strengthen and support the conclusions of the simple analytical model, we 
now present some representative results from our 2.5 D numerical simulation of the general 
two-fluid equations in Cartesian Geometry (Mahajan et al. 2001). For a description of the 
code, the reference (Mahajan et al. 2005) should also be consulted. The simulation system 
is somewhat different because of the existence of an ambient embedding macroscopic field. 
We find that, when such a field is present, the basic qualitatively features of the dynamo 
and reverse dynamo mechanisms do not change much but the algebra is considerably more 
complicated and will be presented in a longer paper later. 

The simulation system contains several effects not included in the analysis; it has, for 
instance, dissipation and heat flux in addition to the vorticity and the Hall terms.The plasma 
is taken to be compressible and embedded in a gravitational field; this provides an extra 
possibility for micro-scale structure creation. Transport coefficients for heat conduction and 
viscosity are taken from (Braginski 1965). 

The simulation presented here deals with the trapping and amplification of a primary 
flow impinging on a single closed-line structure. The choice of initial conditions is guided 
by the observational evidence (Aschwanden et al. 2001; Woo & Habbal & Feldman 2004) of 
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the self-consistent process of acceleration and trapping/heating of plasma particles in the 
hnely structured solar atmosphere. The simulation begins with a weak symmetric up-flow 
(initially Gaussian, where Cgo is an initial sound velocity) with its peak 

located in the central region of a single closed magnetic held structure (location of held 
maximum Bqz = 100 G - upper plot of Fig.l for the vector potential (hux function) dehning 
the 2D arcade). This held was assumed to be initially uniform in time. The magnetic held 
is represented as: 'Q = V x A + with A( 0 ;yly; 0 ); b = B/iloz; bx{t,x,z 7 ^ 0 ) 7 ^ 
0. From numerous runs on the how-held evolution, we have chosen to display the results 
corresponding to the following initial and boundary how parameters: Vornax{,Xoi z = 0 ) = 
Vqz = 2.18 ■ 10^ cm/s; n^max = 10^^ cm“^; T{x,z = 0) = const = Tq = 10 eG. The 
background plasma density is Ubg = 0.2nomax- In simulations n{x,z,t = 0) = n/nomax 
is an exponentially decreasing function of ; 2 . Experience was a guide to for imposing the 
following boundary condition, d^IC^x = ±oo,z,t) = 0 which was used with sufficiently 
high accuracy for all parameters /C(A, T, V, B, n) . The initial velocity held has a pulse-like 
distribution (middle and lower plots of Fig.l) with a time duration to = 100 s. 

It is found that: 

( 1 ) the acceleration is signihcant in the vicinity of the magnetic held-maximum (originally 
present or newly created during the evolution) with strong deformation of held lines and 
energy re-distribution due to magneto-huid coupling and dissipative ehects. 

(2) Initially, a part of the how is trapped in the maximum held localization area, accumu¬ 
lated, cooled and accelerated (plots corresponding to t = 100s in Fig.2). The accelerated 
how reaches speeds greater than lOOkm/s in less than 100 s (in agreement with recent obser¬ 
vations (Schrijver et al. 1999; Seaton et al. 2001; Ryutova & Tarbell 2003) and references 
therein). 

(3) After this stage the how passes through a series of quasi-equilibria. In this relatively 
extended era (~ 1000 s) of stochastic/oscilating acceleration, the intermittent hows contin¬ 
uously acquire energy (see Fig.3 for the how kinetic and magnetic energy maxima and also 
Fig.2 results at t = 1000s). 

(4) The how starts to accelerate again (Fig.3(a-c) for the velocity held evolution). This 
process is completely consistent with the analytical prediction; the acceleration is highest in 
the strong held regions (newly generated, Fig. 2 ). At this moment the accelerated daughter 
hows (macro-scale) are decoupled from the mother how carrying currents and modifying 
the initial arcade held creating new bmax localization areas that span the region between 
< 0.05 and ~ 0.01 Rg from the interaction surface. 

The extensive simulation runs also show that when dissipation is present, the hall term 
(proportion to ao), through the mediation of micro-scale physics, plays a crucial role in 
the acceleration/heating processes. The existence of initial fast acceleration in the region 
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of maximum localization of the original magnetic field, and the creation of new areas of 
macro-scale magnetic held localization (Fig.2, panel for Ay) with simultaneous transfer of 
the magnetic energy (oscillatory, micro-scale) to how kinetic energy (Fig.2, panel for |V| 
and Fig.3 results) are manifestations of the combined ehects of the dynamo and reverse 
dynamo phenomena. The maintenance of quasi-steady hows for rather signihcant period is 
also an ehect of the continuous energy supply from huctuations (due to the dissipative. Hall 
and vorticity ehects). These hows are likely to provide a very important input element for 
understanding the hnely structured atmospheres with their richness of dynamical structures 
as well as for the mechanisms of heating, and possible escape of plasmas. 

Notice, that in the simulation the actual magnetic and generalized helicity densities 
are dynamical parameters. Thus even if they are not in the required range initially, their 
evolution could bring them in the range where they could satisfy conditions needed to effi¬ 
ciently generate hows. The required conditions could be met at several stages. This could, 
perhaps, explain the existence of several phases of acceleration. Dissipation ehects could 
play a fundamental role in setting up these distinct stages; it could, for example, modify 
the generalized vorticity that will hnally lead to a modihcation of held lines and even to the 
creation of micro scales (shocks or fast huctuations). 


3. Conclusions and Acknowledgments 

From an analysis of the two-huid equations, we have extracted, in this paper, the 
’’reverse-dynamo” mechanism — the amplihcation/generation of fast plasma hows in astro- 
physical systems with initial turbulent (micro scale) magnetic helds. This process is simul¬ 
taneous with, and complimentary to the highly explored dynamo mechanism. It is found 
(both analytically and numerically) that the generation of macro-scale hows is an essential 
consequence of the magneto-huid coupling, and is independent of the initial and boundary 
conditions. The generation of macro scale magnetic helds and hows goes hand in hand; 
the greater the macro-scale magnetic held (generated locally) the greater the macro-scale 
velocity held (generated locally). The acceleration due to the reverse dynamo is directly 
proportional to the initial turbulent magnetic energy. When the microscopic magnetic held 
is initially dominant, a major part of its energy transforms to macro-scale how energy; a 
weak macro-scale magnetic held is generated along with. 

The reverse dynamo mechanism, providing an unfailing source for macro-scale plasma 
hows, is likely to be an important mechanism for understanding a host of phenomena in 
astrophysical systems. 
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Fig. 1.— Upper plot: contour plot for the y~ component of vector potential Ay (flux function) 
in the x — z plane for an initial distribution of ambient arcade-like magnetic held. The held 
has a maximum Bmax{xo = 0 , 2 ;o = 0) = 100 G . Middle plot: initial symmetric prohles 
of the radial velocity 14, and density n. The respective maxima (at a:=0) are ~ 2km/s 
and 10^^ . Lower plot corresponds to time evolution of initial how: Vz{t,z = 0) = 

l/o2sin(7rt/to); Vz(t > to) = 0 ; to = 100 s . 
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Fig. 2.— X — z contour plots at 3 time-frames: t = 100 s; 1000 s; 2500 s for the dynamical 
evolution of Ay (first panel from the left), n (second panel), |V| (third panel) and T (last 
panel) for flow - arcade held interaction. The realistic viscosity and heat-hux effects as 
well as the Hall term (ao = 3.3 ■ 10“^°) are included in the simulation. Primary how (type 
displayed in Fig.l) is accelerated as it makes a way through the magnetic held with an 
arcade-like structure (Fig.l). The primary how, locally sub-Alfvenic, is accelerated reaching 
signihcant speeds (> lOOkm/s) in a very short time (< 100 s). Initially the ehect is strong 
in the strong held region (center of the arcade). There is a critical time (< 1000 s) when 
the accelerated how bifurcates in 2; the original arcade held is deformed correspondingly. 
After the bifurcation, strong magnetic held localization areas, carrying currents, are created 
symmetrically about a; = 0. Post-bifurcation daughter hows are localized in the newly 
created magnetic held localization areas. The maximum density of each daughter how is of 
the order of the density of the mother-how. Daughter-hows have distinguishable dimensions 
~ 0.05 At t > 1000 s, the velocities reach ~ 500km/s or even greater (< 800km/s) 
values. The distance from surface where it happens is > 0.01 Rg . In the regions of daughter 
hows localization there is a signihcant cooling while the nearby regions are heated. 
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Fig. 3.— Evolution of maximum values of |V|, \Vp\ = 14 ((a-c)) and 

|b|, l&pl = (&^ + bz ((d-f)) in time. (a),(d) - It is shown that much of the transfer 

from magnetic held energy happens while the hrst and very fast (~ 100 s) acceleration stage; 
(e),(b) - later, the dissipation of perpendicular (towards height) magnetic held huctuations 
lead to the maintenance of the quasi-equilibrium fast perpendicular hows for a period of 
~ 1000s and then the ehective acceleration of how follows; (c),(f) - maximum value of 
magnetic held component along height is not changed and radial component of velocity 
held dissipates ehectively. It should be emphasized that these maximum values of both held 
parameters change the localization dynamically and follow the relationship found analytically 
- fast hows (see Fig.2) are observed in the regions of macro scale magnetic held maximum 
localization (initially given or later generated). 



